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Abstract 

The corrections to the model of nn transitions in the medium and additional baryon- 
number-violating processes are considered. We focus on the time-dependence since it 
is different for the nn transition followed by annihilation (basic process) and processes 
mentioned above. 
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1 Introduction 



In the standard calculations of ab oscillations in the medium [1-3] the interaction of particles 
a and b with the matter is described by the potentials Ua,b- Inif/;, is responsible for loss of 
6-particle intensity. In particular, this model is used for the nn transitions in a medium [4-7]. 
In [8,9] it was shown that one-particle model mentioned above does not describe the total ab 
transition probability as well as the channel corresponding to absorption of the 6-particle. 

In a resent paper [10] the field-theoretical approach to particle oscillations in absorbing 
matter has been proposed. It was studied by the example of nn transitions. In this case 
the main channels are the annihilation and scattering one. The nn transition in the medium 
followed by annihilation 

(n — medium) — >• (n — medium) —>■ M, (1) 

(M are the annihilation mesons) is shown in Fig. la. The nn transition with n in the final 
state is shown in Fig. lb. The block T^^ involves all the n-medium interactions followed by 
annihilation including the antineutron rescattering in the initial state; S — 1 + iT. Similarly, 

involves scattering proper and annihilation loops. Due to this the antineutron propagator 
is bare. This model is considered in more detail below. 

Let (Ta and cr^ are the cross sections of free-space nN annihilation and nN scattering, 
respectively. At the low energies (Jq > 2.5(75 [11,12]. Hence the annihilation channel (Fig. 
la) is dominant: in the first stage of n-medium interaction the annihilation occurs. At the 
nuclear densities only the diagram la is essential because the antineutron annihilates in a time 
Ta ~ 1/r, where F is the annihilation width of n in the medium. 

In [10] the diagrams 1 have been calculated. The corrections to this model and other 
baryon-number-violating (AB = 2) processes [13] (see Figs. 2 and 3) can essentially change 
the result. Besides, it turns out that the time-dependence of the processes shown in Figs. 1 
and 2 is different. This is non-trivial circumstance and we focus on this point. The heart of 
the problem is as follows. The processes shown in Fig. 2 are described by the exponential 
decay law. The diagram la contains the infrared divergence conditioned by zero momentum 
transfer in the nn transition vertex. This is unremovable perculiarity. The fact that amplitude 
is singular means that the standard (S-matrix approach is inapplicable. As a consequence the 
other surprises take place as well, in particular, the different functional structure of the result 
and non-exponential behavior. These problems are studied below. 
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2 Calculation 



First of all we consider the incolierent contribution of the diagrams 2. In Fig. 2a a meson is 
radiated before the nn transition. The interaction Hamiltonian has the form 

Hl^g'^t^'^n + nnn + n, (2) 
Hnn = e*^*n + H.C. (3) 

Here Ti-nn is the Hamiltonian of the free-space nn transition [6,10], e is the off-diagonal mass 
(transition mass); e = 1/Tnn, where Tnn is the free-space nn oscillation time; Tl is the Hamil- 
tonian of n-medium interaction taken in the general form. In the following the background neu- 
tron potential is omitted; m„ = rUn = m. The neutron wave function is np{x) = il^^^"^ exp(— ipx), 
were p = {po, p) and po = m + p^/2m. 
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Figure 1: a. nn transition in the medium followed by annihilation, b nn transition in the 
medium with n in the final state. 

For the process amplitude one obtains 

M2a = 9GeGM^"-'\ (4) 

Q _ ]. /r^\ 

Po — qo — m — {p — q)^/2m + iO' 
where q is the 4-momenta of meson radiated, M^^^ is the amplitude of antineutron annihilation 
in the medium in the (n) mesons. It is defined in the usual manner: 

iTf- =<(n)0| Texp(-i J dxH{x)) - 1 |Onp_g>= iV(27r) - pi)M('*), (6) 
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5* = 1 + iT. Here | Oup^g > is the state of the medium containing the n with the 4-momentum 
p — q, < (n) I denotes the final state of the annihilation mesons, N includes the normalization 
factors of the wave functions. Since M^^^ involves all the n-medium interactions followed by an- 
nihilation including the antineutron rescattering in the initial state, the antineutron propagator 
G is bare; the n self-energy E = 0. 
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Figure 2: Corrections to the model shown in Fig. la (a and d) and additional baryon-number- 
violating processes (b and c). 

Let and F^") be the widths corresponding to the Fig. 2a and annihilation width of n in 
the (n) mesons, respectively; F = J2{n) T^"^- Taking into account that F^") is a smooth function 
of y/s and summing over (n), it is easy to get the estimation: 

F2„ 5 • 10-V^r ^ ^F. (7) 

The time-dependence is determined by the exponential decay law: 

W^2a(t) = l-e-^^"*~F2„t (8) 

If g — > 0, the amphtude increases since G ^ Gg, 

Po — m — p^/2m 

(The limiting transition g — > for the diagram 2a is an imaginary procedure because in the 
vertex n —>■ n$ the real meson is escaped and so go ^ The fact that the amplitude 

increases is essential for us because for Fig. la g = 0. 

Consider now the baryon-number-violating decay n n$ [13] shown in Fig. 2b. It leads to 
the same final state, as the processes depicted in Figs, la and 2a. As with Tinn, the Hamiltonian 
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of the decay n —>■ n$ is taken in the scalar form e$^^$^„+H.c. Then 

ni = e^^+<^>^n + n, (10) 

e$ is dimensionless. The process amphtude is 

M26 = e^GM^'^-i). (11) 

The antineutron propagator G is given by (5). Denoting | q |= for the decay width one 
obtains 

Qo — q'^ + "^1- As with Fig. 2a, the time-dependence is determined by the exponential decay 
law W2b{t) ~ r2bt. 

In Fig. 2c the baryon-number-violating conversion n — > A in the medium [13] is shown. 
The amplitude is 

M2e = eA — 7Mi"\ (13) 

Po-mA- p^/2mA + ^0 

where sa corresponds to the vertex n — > A, is the amplitude of A annihilation in the 

medium in the (n) mesons. Contrary to Fig. la, there is no infrared singularity because 
rriA 7^ m. We also note that this process cannot produce interference, since there is X-meson 
in the final state. 

Let us suppose that there is residual n-medium interaction V (scalar field) which cannot 
be involved in H, i.e. the block T^- (see Fig. 2d). In this case it is included in the in medium 
antineutron Green function 

G„ = G. + G,VG, + ... = py5^ = 4 (14) 

The process amphtude is — eG^M'^'^^ . The process width V2d is found to be 

= ^T. (15) 

In fact, the correct definition of antineutron annihilation amplitude (see (6)) excludes the 
diagrams like 2d. It was adduced only for the purpose of illustration since the models with the 
dressed antineutron propagator were at one time discussed. It also demonstrates the solution 
stability to the possible perturbations, which seems important. 

In Figs. 2c and 2d the vertexes n ^ K and n ^ n are depicted by 2-tail diagram. Neverthe- 
less, the amplitudes M2c and M2d are non-singular because itla 7^ ttl (Fig. 2c) and = m + V 
(Fig. 2d). One can say that effective momentum transfer qo = niA — m and = V takes place. 
The nonzero momentum transfer is the common property of the diagrams shown in Fig. 2. 



5 



n 



n 



n 



n 



n 



n 



-u. 



Figure 3: Free-space nn transition (a) and medium correction to the nn transition vertex (b). 
The hmit obtained for the free-space nn oscillation time is [10,14] 

Tnn > 10'^ yr. (16) 

In [10] and calculations made above the free-space nn transition operator has been used. This 
is impulse approximation which is employed for nuclear /3 decay, for instance. The simplest 
medium correction to the vertex (or off-diagonal mass, or transition mass) is shown in Fig. 3b. 
In this event the replacement should be made: 

e^e^ = e(l + Ae), (17) 

= esb/e, where 63;, is the correction to e produced by the diagram 3b. Then the limit becomes 

r„^>(l + Ae)10i^yr. (18) 

Obviously, the Ae cannot change the order of magnitude of Tnn since the n — > ri operator 
is essentially zero-range one. The free-space nn transition comes from the exchange of Higgs 
bosons with the mass mn > 10^ GeV [5]. Since niH ^ rnw {iriw is the mass of 1^-boson), the 
renormalization effects should not exceed those characteristic of nuclear /3 decay which is less 
than 0.25 [15]. So the medium corrections to the vertex are inessential for us. 

We now turn to a discussion of interference between different processes. The diagram 2c 
cannot produce interference, since it contains X-meson in the final state. Figure 2d is included 
in diagram la. As mentioned above, it was adduced only for illustration of the effect on our 
results of possible fluctuations. For the rest of the diagrams the significant interferences are 
unlikely because the final state in nN annihilation are very complicated configurations and 
persistent phase relations between different amplitudes cannot be expected. This qualitative 
picture is confirmed by our calculations [16] for p-nuclear annihilation. Moreover, in [16] the 
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interference between the amplitudes which give the comparable contribution takes place. As 
shown later, the contribution of diagrams 2 is negligible. 

The corrections given above should be compared with the results obtained for the basic 
model. We adduce the main results corresponding to diagram la. (This problem has been con- 
sidered in [10]. Nevertheless, we should make additional comments since the time-dependence 
of the processes shown in Figs. 1 and 2 is different (see below) what is very unusual.) The 
interaction Hamiltonian is taken in the general form: 

Hi = Hnn + n. (19) 

The notations are the same as in (2). Formally, in the lowest order in T-Lnn the process amplitude 
Ml is uniquely determined by the Hamiltonian (19): 

Ml = eGsM^""^ (20) 

(see Fig. la), where the antineutron propagator Gs is given by (9). So Mi ~ 1/0. This is 
infrared divergence conditioned by zero momentum transfer in the nn transition vertex. For 
solving the problem the approach with a finite time interval (FTA) [17] is used. The problem 
is formulated on the interval (t, 0) which corresponds to the concrete conditions of the process 
under study. (Strictly speaking diagrams 2 should be calculated on the interval (t, 0) as well. 
However, for non-singular diagrams the FTA converts to ^'-matrix theory [10].) If 

n > 1 (21) 

then the process (1) probabihty Wi{t) is 

Wiit) ^ Wf{t) = eH\ (22) 

where Wf{t) is the free-space nn transition probability. Owing to the strong annihilation 
channel, Wi{t) ^ Wf{t). The antineutron annihilates in a time Tq ~ l/F. So condition (21) 
implies that the annihilation can be considered instantaneous: Ta ~ l/F <^ 

Equation (22) corresponds to the limiting case (21). The opposite hmiting case F — > 
means that there is no annihilation and so l^i = 0. This result follows from the equations of 
motion. In the intermediate range Fi ~ 1 the distribution Wi{t) is more complicated (see Eq. 
(64) and Sect. 8 of Ref. [10]). It can be expected that Wi{t) < Wf{t), although no detailed 
calculations is possible. 

The physical models used for diagrams 1 and 2 are identical. The sole (but fundamental) 
distinction of Fig. 1 with respect to the diagrams 2 is the zero momentum transfer in the first 
vertex. Also we would like to emphasize that result (22) is true for any neutron wave function 
since the zero momentum transfer in the nn transition vertex takes place in any case. 
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3 Time-dependence and contribution of corrections 



Comparing (22) and (8) we have 



W2a{t) 

where we have put m$ = mjr. Consequently, if 



r 



(23) 



mlt/T » 1, 



(24) 



and rt ^ 1 (see (21)) then the contribution of diagram 2a is neghgible. (For the nn transition 
in nuclei conditions (21) and (24) are fulfilled since in this case F ~ 100 MeV and t = Tq = 1.3 
yr, where Tq is the observation time in proton-decay type experiment [18].) In fact, it is suffice 
to hold condition (21) only because it is more strong. 

Similarly, with the help of (22) and (15) we conclude: if V'^t/T ^ 1 and Ft ^ 1, the 
contribution of diagram 2d is negligible. So ^ Wi and ^ Wi because ~ t 
and W2d ~ t, whereas Wi ~ t^. The quadratic timc-dcpcndcncc of the Wi is due to the zero 
momentum transfer in the nn transition vertex. The heart of the problem is as follows. 

For the processes with q = the ^'-matrix problem formulation (oo, — oo) is physically 
incorrect [10]. This is apparent even from the limiting case 7i = 0: if Tij — H-nn (see (19)), 
the solution is periodic. It is obtained by means of non-stationary equations of motion and 
not S'-matrix theory. To reproduce the limiting case 7i — > 0, i.e. the periodic solution, we 
have to use the FTA. Formally, the <5'-matrix approach is inapplicable because the perculiarity 
Ml ~ 1/0 is unremovable one. The FTA is infrared-free. 

If the problem is formulated on the interval (t, 0), the decay width F cannot be introduced 
since F = S/^i I -S'/j(oo, — oo) |^ /Tq, Tq — > oo. This means that the standard calculation 
scheme should be completely revised. (We would hke to emphasize this fact.) The direct 
calculation by means of evolution operator gives the distribution (22). 

The more physical explanation of the t^-dependence is as follows. In the Hamiltonian (2) 
corresponding to Fig. 2a we put li — Hnn — 0. Then the virtual decay n — > n$ takes place. 
The first vertex of the diagram 2a dictates the exponential decay law of the overall process 
shown in Fig. 2a. Similarly, in the Hamiltonian (19) corresponding to Fig. la, we put H — 0. 
Then the free-space nn transition takes place which is quadratic in time: Wf{t) = eH'^. The first 
vertex determines the time-dependence of the whole process at least for small F. (See, however, 
text below (22).) We also recall that even for proton decay the possibifity of non-exponential 
behavior is realistic [19-22]. 

In Fig. 2c the transition n ^ A is depicted by 2-tail vertex. Nevertheless, the diagram 2c 
is described by the exponential decay law because in the above-mentioned vertex the effective 
momentum transfer takes place. As a result, the structure of the amplitude M2c is the same as 
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that of the and M2b- Indeed, let us compare (4) and (13). Put p = q = for simpUcity. 
Equation (5) becomes G = l/{po — qo — m + iO). On the other hand, if Qq = m\ — m, this 
expression is the propagator of (13). Since the structure of propagators and amplitudes 
and is identical, the M2c corresponds to the effective decay. Due to the effective momentum 
transfer q — {m^^ — m, 0) in the first vertex, the diagram 2c is described by the exponential 
decay law as well. 

Ratio (23) does not depend on e, whereas W2b/Wi contains two unknown parameters e and e$ 
and so the contribution of diagrams la and 2b cannot be compared numerically. Nevertheless, 
it is safe to argue that for sufficiently large times W2b <^ Wi because Wi{t) ~ i^, whereas 
W2b{t) ~ t. The same is true for Fig. 2c: W2c < Wi. 

Formally, the different time-dependence of the processes shown in Figs. 1 and 2 is due to 
^'-dependence of amplitudes. If q decreases, the amplitudes M2a — M2d increase; in the limit 
g — > they are singular. The point q — corresponds to realistic process shown in Fig. la. The 
i^-dependence of this process is the consequence of the zero momentum transfer. Due to this 
for sufficiently large times {t ^ l/T for Fig. 2a, for example) the contribution of diagrams 2 is 
neghgible. We also recall that the diagrams 2c and 2d give no contribution to the amphtude 
of the process under study. They have been considered for illustration of the role of additional 
baryon-number- violating processes (Fig. 2c) and possible fluctuations (Fig. 2d). 

4 Conclusion 

The diagrams shown in Fig. 2 correspond to the decays (Figs. 2a and 2b) or decay-type pro- 
cesses (Figs. 2c and 2d). They are described by the exponential decay law. The corresponding 
decay widths have been calculated. The probability of the process shown in Fig. la is quadratic 
in time: Wi ~ t'^ (see, however, text below (22)). This is due to the zero momentum transfer 
in the nn transition vertex. The point = is the peculiar point of the (S'-matrix amplitude. If 
the amplitude is singular, the surprises can be expected. From this standpoint a departure from 
the exponential decay law comes as no surprise to us. It seems natural that for non-singular and 
singular diagrams the functional structure of the results is different, including t-dependence. 
The opposite situation would be strange. 

To summarize, the diagrams la and 2 are fundamentally different. The time-dependence 
of the corresponding processes is different as well. Because of this for sufficiently large times 
the contribution of the diagrams 2 is neghgible. The distribution (22) and hmit (16) can be 
essentially changed only if there is no diagram la, which is out of question. For the K'^K^ 
oscillations in the medium the similar consideration can be done. 
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